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Method to estimate the parameters of a propagation channel 

BACKGROUND OF THE INVENTION 
1. Field of the Invention 

The invention relates to a method for estimating the impulse response of 
5 a propagation channel, especially its parameters such as the delay values, the 
directions of arrival as well as the values of the complex amplitudes associated with 
these parameters, with an a priori knowledge of the signal. 

It is applicable, for example, to the estimation of the parameters of a 
finite impulse response filter which can be written, though not necessarily so, in 
10 specular form, namely a filter that can be written in the form of weighted Diracs. 

The invention can be used to estimate propagation channels in the field 
of radiocommunications but also, generally, it can be applied to any signal filtered by 
a finite impulse response linear filter. 

In a transmission system, especially one using radio waves, a transmitter 
15 sends out a signal in a transmission channel to a receiver. The signal that is sent 
undergoes amplitude and phase fluctuations in the transmission channel. The signal 
received by the receiver consists of copies of the transmitted signal that are 
temporally shifted and modified. The fluctuations of the signal and the shifts give 
rise to a phenomenon known to those skilled in the art as inter-symbol interference. 
20 The interference arises especially out of the law of modulation used for transmission 
and also from the multi-path propagation in the channel. 

The received signal generally arises out of a large number of reflections 
in the channel. Hie different paths taken by the sent signal thus lead to different 
delays in the receiver. The impulse response of the channel represents the totality of 
25 the fluctuations to which the sent signal is subjected. 

The estimation of the propagation channel in a radiocommunications 
system is useful in several respects, some of which are indicated here below as 
examples. 

• The demodulators generally require knowledge of the channels in order to 
30 remedy the harmful effects that they have caused, 

• a second point of interest is that of urban or extra-urban localization, for example 
the principle of emergency localization using the number "91 1 " in the United States 

• Finally, the knowledge of the propagation channels can also serve for the use of 
smart antennas at reception as well as transmission. 

35 2. Description of the Prior Art 
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There are various techniques known in the prior art for estimating 
propagation channels and their parameters. 

For example the document by R. Rick and L. Milsteil, "Performance 
acquisition in mobile ds-cdma systems", IEEE Trans on Communications, Vol: 45 
5 (No; 1 l):pp: 1466-1476, November 1997, proposes a search for propagation delays 
by using a bank of non-coherent detectors. The results are proposed for multi-path 
channels in the presence of Doppler phenomena and inter-cell and intra-cell 
interference. 

The document by R. Rick et L- Milsteil, "Optimal decision strategies for 
10 acquisition of spread-spectrum signals in frequency selective fading channels" in 
IEEE Trans* on Communications, Vol: 46 (No: 5); pp: 686-694, May 1998, discloses 
an optimal decision rule based on the outputs of the correlators proposed in the 
document referred to here above. A single-user technique of this kind is limited by 
interference in the case of multiple-users. 
15 There also exist known ways of using rectangular shaping filters, for 

example by the method described in one of the following documents: 

• E. Strom, S. Parkvall, S. Miller, and B. Ottersen, "Propagation delay estimation 
in asynchronous direct-sequence code-division multiple access systems", IEEE 
Trans on Communications, Vol: 44:pp: 84-93, January 1996 

20 • S. Parkvall, "Near-Far Resistant DS-CDMA Systems : Parameter estimation and 
Data Detection", PhD thesis, Royal Institute of Technology Stockholm, Sweden, 
1996. 

• S, E> Bensley and B. Aazhang, "Maximum likelihood estimation of a single user's 
delay for code division multiple access communication systems", Conf. Information 

25 Sciences and Systems, 1994. 

In the case of shaping filters with a duration greater than a chip time, 
these different methods are no longer suitable. 

Algorithms for the combined estimation of angles of arrival and of the 
differential delay times, on known and received signals, based on sub-space 
30 techniques have been proposed, for example in the document by P. Gounon, 
"Analyse spatio-temporelle haute resolution k Taide d'une antenne active", (High 
Resolution Space-time Analysis Using an Active Antenna ) Traitement du Signai 
(Signal Processing), Vol. 1 1 (No. 5), pp. 351-360, 1994. 

The document by A. J Van der Veen, M* C. Vanderveen, et A. J. Paulraj, 
35 "Joint angle and delay estimation using shift-mvariance properties", IEEE Sig. Proc 
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Letters, VoL4 (No.5): pp. 142-145, 1997, discloses methods for the estimation of the 
physical parameters of propagation by means of methods based on sub-spaces. 

However, such methods suffer from a deterioration of performance 
characteristics once the impulse responses of the propagation channels are correlated. 
5 This situation occurs especially when the complex amplitudes do not vary with 
sufficient speed on the covariancc matrix of the impulse responses estimated in terms 
of the least-error squares by means of the signal transmitted. 

Figure 1 shows the different techniques of maximum likelihood. 
A maximum likelihood method has been proposed, for example, in one 
1 0 of the following references: 

• J. Grouffaud, M Identification spatio-temporelle de canaux de 
propagation & trajets multiples", (Space-time Identification of Multi-Path 
Propagation Channels), PhD thesis, Ecole Normale Superieure de Cachan, June 
1997. 

15 • M. Wax and A. Leshem, "Joint estimation of delays and directions of 

arrival of multiple reflections of a known signal.", IEEE Trans, on Signal 
Processing, Vol: 45(No: 10):pp; 2477-248, October 1997. 
but it does not deal with the MIMO (Multiple Input Multiple Output) context. 

The document by R Graffoultere, "Methodes actives spatio-temporelles 

20 large bande" (Active Wideband Space-Time Methods), published in « Techniques et 
performances. Applications En Sonar » (Techniques and Performance. Sonar 
Applications), PhD thesis, INPG, March (997, also discloses u method of estimation 
based on maximum likelihood but the studies on performance deal only with the case 
of a single source or of several distinctly separate sources. A similar study is 

25 disclosed in the document N. Bertaux , "Contribution a 1'utilisation des methodes du 
Maximum de Vraisemblance en traitement radar actif 1 (Contribution to the Use of 
Maximum Likelihood Methods in Active Radar Processing), PhD thesis, Ecole 
Normale Superieure de Cachan, January 200G, lor active radar applications in the 
case of single sources. 

30 SUMMARY OF THE INVENTION 

The present invention relates to a method that can be used especially to 
estimate the parameters of the propagation channel by working on the correlated 
signals, in selecting a certain number of samples and in searching for the values of 
the delay parameters and/or directions of arrival, for example, which would enable 

35 the most efficient reconstruction of the signal received. 
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The invention also relates to a method that integrates pulse compression 
techniques such as pre-processing. 

The invention relates to a method for the estimation of one or more 
parameters of a propagation channel with a priori knowledge of the signal in a 
5 system comprising one or more sensors. 

The method of the invention comprises the following steps: 

• correlating the signal or signals x(t) received by the sensor or sensors with 
a known signal c(t), 

• sampling said signal after correlation at a sampling period Te and selecting 
10 a number of samples per concatenation, 

• determining at leas tone parameter of the propagation channel such as t 
and/or 0 which enables the most efficient reconstruction of the signals received by 
using a maximum likelihood method. 

According to one embodiment, the characteristics of the system of 
15 sensors are known and comprise for example: 

• a correlation step with a known signal c(t) equal to 1 , 

• the signals received on the antenna being expressed in the form X=S(0/r)h+B 7 
and 

• the estimates of the parameters as T and 0 being expressed in the following form: 

20 

d t ? = argminjjn£(0,T)Xf 
- arsminfX^C^rJX} 

25 where n^s is the projector orthogonal to the image generated by the 

column vectors of S(t, 0). 

The characteristics of the system of sensors are for example the response 
to the antenna. 

It may comprise a srep for determining the complex amplitudes h of the 
30 impulse response of the propagation channel from the estimates of the estimated 
parameters as t and 0. 

According to another embodiment, the characteristics of the system of 
sensors are not known, and the method comprises for example: 

• a step for the correlation of the signals received by the network of sensors 
35 with a known signal c(t) equal to 1 , 



♦ a step in which the signals received are expressed in a concatenated form 
Y 2= \|/(t)cw-N where y(T) is equal to the convoluted product of the unit matrix In 
and the matrix 5(T)^ , (r]),»^(Tpp,.--,j(r^)J J and a contains the responses of the 

paths of the different users, 

• a step for the estimation of the delay vector t from 

r = argmm |jni(r)Y|| 
a= argmin/r 

fytni(r)Y) 

where TT\j7 is the projector orthogonal to the image generated by the iine vectors of 

The step of correlation of the signals is performed with a signal c(t) 
different from L The characteristics of the system of sensors are known and this 
correlation step comprises a step for the estimation of the parameters t and 9 from: 



where 



0 7 r m argmin^R^nWrjX" 
n£ - 1 - *(*, r) (*t(fi, T )B£ r)) ~ l r)B£ 1 



The invention also relates to a device for estimating one or more 
parameters of a propagation channel with a priori knowledge of the signal in a 
system comprising one or more sensors. 

The device comprises at least; 

• a device adapted to the correlation of the signal received by the sensor or 
sensors s(t) with a known signal c(t), 

» a device adapted to the selection of a number of samples of the signal 
obtained after the correlation step, and 

• a device adapted to the determining of the parameters of the channel by a 
maximum likelihood method. 

The method can be applied, for example, in applications of the MIMO 
(Multiple Input Multiple output) or SIMO (Single Input Single Output) type. 
The invention in particular has the following advantages: 
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• in the cooperative MJMO context, the estimator, within the framework of a 
known antenna, makes use of a parametric model of the signals received on the 
antenna, unlike the deterministic maximum likelihood method which estimates each 
sample of the signal without any a priori assumptions, 
5 • in limiting the field of analysis in certain cases of application, it leads to a 

reduction of the size of the data processed, and thus results in faster processing. 
BRIEF DESCRIPTION OF THE DRAWINGS 

Other features of the invention shall appear from the following 
description which provides examples of implementation given by way of an 
10 illustration, with reference to the appended drawings, of which; 

• Figure 1 shows different approaches of estimation of parameters using the 
Maximum likelihood method 

• Figure 2 gives a schematic view of a first alternative embodiment of the 
method according to the invention in which the characteristics of the reception 

15 antennas are known* 

• Figure 3 is a drawing of a second alternative embodiment in which the 
characteristics of the reception antennas are not known, 

• Figure 4 gives a schematic view of another variant of the method 
comprising a preliminary step that can be applied to the methods described in figures 

20 2 and 3 here above. 

MORE DETAILED DESCRIPTION 

In order to provide a clearer understanding of the object of invention, the 
description is given by way of an illustration that in no way restricts the scope of the 
invention in the field of mobile radio communications, in a cooperative context in 
25 which the receiver has certain information available to it, for example the 
knowledge, to a certain extent, of the signal sent (a priori knowledge of the signal) 

The receiver is constituted for example by N identical and 
omnidirectional sensors. The signals s(t) received by each of the sensors are 
regularly sampled at the rate T e . This example does not prevent the application of the 
30 method of the invention to other types of receivers. 

Before giving a detailed description of the different alternative forms of 
implementation of the method, we shall recall a few points concerning the modelling 
of the signals and shall state a few assumptions enabling the implementation of the 
method. 

35 Modelling of signals received on the reception system 



When there are U users (sources), the U .signals s u (t) with 1 ^ " £ U are 
propagated through different channels before being received by a set of sensors 
indexed k. The impulse response of the channel between the user u and the sensor k 
is referenced h\(t, t) where t is the instant of arrival of the signal and x is its delay. 

If we call h\(U x) the vector with the dimension {Nx\) whose 

components are the values h\(U T), with I < k < K and K the number of sensors, the 

models of the received signals can be written in the form (1): 

u 

x(t)=Xs w (T) *h U (t > t)+b(t) where 

tfa) 

b(t) is the noise vector with the dimension (Nxl). 
Specular Assumption 
Time-related specular assumption 

The specular model is a widely used model of channels for which it is 
considered that the paths are perceived by the antenna as being temporally discrete 
paths. 

The transfer function of the channel therefore becomes (2): 
m.r> X h u p (T p «(t))5(t^ p "(t)) 

with 

• Pu the number of propagation paths of the user u , 

* h p (Tp u (t)) lhe {Nx\) vector containing the product of the responses of the 

antenna to the path p multiplied by the complex amplitude of the path. 
The received signal can be written in the case of the specular channej (3) 

Calibrated antenna 

If the antenna is calibrated, assuming that the signals reaching the 
antenna have a negligible angular dispersion and the assumption that the signals are 
narrowband signals is met, the received signal is modelled by (4): 

*(') - EE^W^' " r " (0) + m 

where 

0" is the direction of arrival of the path p of the user u, 
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U 

• p ^ T p U ^ * s ^ e cor npIex gain of the path p of the user u received on the first 
sensor. 

The previously made assumption (concerning the specularity of the 

direction of arrival) is not restrictive in the case of a base station in height. If the 

5 specularity is not sufficient, the model may be modified by adding a lobe width in 

the direction of arrival which will be estimated by the method proposed in the 

invention by adding the parameters of this lobe to the parameters to be estimated. 

The different variants of the method according to the invention use the 

maximum likelihood method, for which various approaches to the estimation of 

10 parameters have been given in figure 1 . 

First alternative embodiment: calibrated antenna 

A first alternative embodiment of the method according to the invention 

uses a maximum likelihood method adapted lo the signals sent, when these signals 

sent are partly known, and the antenna response for all the directions of arrival is 

15 known. 

Figure 2 shows a block diagram where the processing block I, 
comprising for example a microprocessor, adapted to implementing the method 
according to the invention, receives firstly the different signals received by the 
antennas 2i as well as a known signal corresponding to the signal sent. 

20 At output of the processing block, the method according to the invention 

has at least estimates of the direction of arrival 8 and delay vectors x which, for 
example, are used in order to deduce therefrom the impulse response of the channel 
h(t), for a given user or for ail the users. 

The signals received on the antenna array are delayed, expanded and, 

25 possibly, frequency-transposed versions of the signals sent by the users present in a 
cell (or propagation channel). These signals are therefore expressed as a linear 
combination of a set of known basic functions whose parameters are to be extracted. 

The method is presented in the general context of parametric estimation 
of the angles 9 and of the delays T for any antenna whatsoever constituted by N 

30 identical and omnidirectional sensors, The signals are sampled regularly at the rate 
T e . It is assumed that the number of paths are known and that, if this is the case, the 
number is determined by using a method of detection of the number of sources by 
which is possible to obtain the number of paths, described for example in the 
following references. 



• J. J, Fuchs, "Estimation of the number of signals in the presence of unknown 
correlated sensor noise", IEEE Trans on Signal Processing, Vol 40 (N° : 5) pages 
351-360, 1994; 

• J.Gouffraud* P.Larzabal and H.Clergeot, "Methodes h Haute resolution pour le 
traitement d'antenne et r analyse spectrale" (High-resolution methods for antenna 
processing and spectral analysis) in « Selection d'ordre de modele », Editions 
Hermes, February 1998, Chapter 6; 

• H. Akaike, "A new look at the statistical model identification IEEE Trans. On 
Automatic Control, Vol 19 (No. 6), December 1974, 

• M. Wax and T. Kailath, "Detection of signals by information theoric criteria" de 
M.Wax et T.Kailath, IEEE Trans.on ASSP, Vol 33(N Q 2), pages 387-392, April 
1985. 

The steps implemented in the algorithm according to the invention 
remain valid when the number of sources is overestimated. 
Modelling 

The signal x(t) (equation (4) received by the array of sensors 2i is 
sampled at the rate T €y during a period t=[T e M e T e ] with Ne being the number of 
samples or number of observations during which the complex amplitudes of the 
impulse response for a user u and a path p, h^ix^it)) are considered to be 

constant. 

These N„ observations with the dimension Nx] are concatenated in the 

t T T 

vector X with the dimension N e Nx\:X=[ x(T e ) V-%X(N e T e ) J (5) where the 
exponent T designates the transpose. 

We define the vector B with the dimension N e Nx\ containing the 
concatenation of the samples of the noise during the same period by (6): 
BKb(T e ) T "%b(N e Te) T ] T 

Finally, s w (T p u ) is taken to denote the vector with the dimension N e x 1 
corresponding to the signal sent by the user a, delayed by (x p u and sampled: 

s M «VW s U (Tc-(V)^-^ U ( N eT e -(V)] T (?) 

The vector of the concatenations of the signals received can therefore be 

written as follows (8): 

u p m 

X=Z I h u ^x u p )(s u {x u p )®ci(ep)^B 
where ® is Kronecker's product 
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Let S(8,T) be the matrix N e NxP (p=^p u ) of the delay-direction vectors 
defined as follows (9) : 

m r) -[[.-«)« ■<*] . • • • ■ [■" W • <*>] ■ • • ■ ■ » • )l] 



where: 

- 8 is the F-sized real vector of the directions of arrival. 

- x is the P-sized real vector of the delays, (P being the number of 
10 paths) 

The equation of the signals received on the antenna can therefore be 
written as; X=S(0,T)h-riB (10) where h is the Pxl complex vector of the amplitudes. 

This formulation has the value of expressing the observations us a linear 
model that is a function of the complex amplitudes h and of the noise. 
15 It must be clearly noted that this is a MIMO (multiple input multiple 

output) modelling, grouping together all the signals of all the users. 

It is not necessary to know all the signals sent. For example, in mobile 
communications, sequences of known signals are sent periodically, resulting in a 
known signal. In the radar or in radio altimetry, for example, the transmitted pulse is 
10 known. 

Log-likelihood 

The complex amplitudes h are considered to be unknown but 
deterministic. In order to simplify the explanation of the invention without, however, 
limiting it, the signals s ll (t) forming S(6,t) are considered to be perfectly known, 
S apart from the parameters 0/c. This leads to an observation model (10) where only 
the additive noise vector is random, with a Gaussian probability density, represented 

by (II): 

P {X|«r 2 , h,' 0, r } = ff „.„* 3/v .„ exp (- £ (X - S(0, r)h>< {X - S(0, r)h)) 



The log»iikelihood of the observations is therefore (12): 

L (Xlo 2 , h, *, r) * -KN Jog(jr^) - «X - S(9, r)hf 
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Estimator of the parameters tf 2 , h, 0, % 

Through the separability of the parameters a 2 , h, % x, the expression of the 
log-likelihood (12) authorizes the analytical determining of the vector of the complex 
amplitudes h and of the power of the noise a 2 as a function of 0 and T. These 
5 estimates are both obtained by cancelling the derivatives of the log-likelihood 
respectively with respect to a 1 and with respect to h. Their expressions are 
respectively (13) and (14); 



10 and 



fi § 

jft nil 



By substituting the values of h and of a 2 in the log-likelihood (12), the 
1 5 expression gets simplified and the estimates of 8 and x are given by ( 1 5) 



0, t - argmin |Jx - S(0, r) (S*(«, r)S(0, r)) ~ ! S»(0, r)X| 



20 In letting Tlj denote the projector orthogonal to the image generated by 

the column vectors of S(8, T) ( 16) 

n x (ff) r) , 1 3{9t r) (S»(*. r)S(«, r))" 1 S« (ff. r) 

25 

The estimates are the solutions of (17) 

S,r = argmin||n^(e,r)Xj| 2 

30 

These estimated values are then used respectively in the expressions 
(13) and (14) to determine the power of the noise and the impulse response of the 
propagation channel. 
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The criterion introduced in (17) is referenced hereinafter in the text by: 



Implementation of the method according to the invention in the case of a non- 
5 calibrated antenna 

Figure 3 gives the schematic view of an exemplary implementation of the 
method in the case where the characteristics of the receiver or the antenna are not 
known. The antenna is formed by several sensors, whose response is not known, and 
is often designated by the expression "unstructured antennas". 
10 The antennas 2i receive the signals that are transmitted to a processing 

block 14 adapted to implementing the steps of the method, this block itself receiving 
a known signal 

The processing block, which is equipped for example with a 
microprocessor programmed accordingly, delivers at least the estimated value of the 
15 delay vector t and the estimated value of the impulse response of the propagation 
channel* 

Modelling 

The received signal is expressed in the form of the equation (4) and is 
sampled at the rateTe on a period of time t = [Te, NeTe] in which the complex 
20 amplitudes of the impulse response of the propagation channel hp U (tp U (0) may be 
considered to be constant or substantially constant* These Nc samples with a 
dimension Nxl are concatenated in a vector Y with a dimension NeNxl such as 
(18): 



25 



30 



Y = 



JD(7i)„.. t to(N*T<\ t ..,XN(T<>) f ...,XN(N<T:) 



Minor t 



scn,suf N 



Where xi(nTe) is the sample n of the sensor i. 

Similarly, the NeNxl sized vector N contains the concatenation of the noise samples 
(19). 



- — v — ~ — 



Where bi(nTe) is the noise sample n of the sensor with an index i* 
The signal Y can be written (20) 
Y = \|/(TT)a+N 
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where the matrix (N^xNP) >j/ is buik from the different values s"(T p u ) as follows 
(21): 

y(x) = In<8>S(t) 
where 5( T ) is the matrix (N„xP) such as (22): 

a contains the responses of the sensors to the different paths of the different 
users (23) : 

if 



CL = 



sensor I 



with: 



The complex amplitudes are always assumed to be constant during the phase 
of estimation of the parameters. 

The modelling of the received signals given by the equation (20) provides us 
with an affine model . As above, this expression takes account of the signals of all 
the users. 

Log-likelihood 

The complex amplitudes a are considered in this case as being unknown but 
deterministic , The signals that come into play in ¥(t) are assumed to be known apart 
from the variables x . This leads to a model where only the noise N is random, with 
Gaussian components and with a probability density (25) : 

P { Y|<r*, «, r } = exp ((Y - *(r)a)'(Y - *(r)«))) 

The log-likelihood is then given by (26) 

Z(Y|<r>, a, t) - -N e N tatfw 3 ) - ~ | Y - *(r«)!| 2 

Estimator of the parameters cr 2 , a, x 

As in the case of the maximum likelihood adapted to the signals received 
with a known antenna as explained here above, the analytical expressions of the 
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complex amplitudes a as well as the power of the noisea 2 may be expressed as a 
function of the set of parameters of the delay vectors t. 

This can be obtained for example by respectively cancelling the derivatives of 
the log-likelihood with respect respectively to a and to cr 2 > 

We obtain the expression (27) 

and (28) 

« 



a ~ (» t (r)«<r)p ¥ f (r)Y = **(r)Y 

By replacing a and <fi by their estimates, the log-likelihood gets simplified 
and the estimator of X is (29) 

r = argnua IJy - #(r) (•'M#(r))~ l * r {r)Yj| 2 

If n^(T) denotes the projector that is orthogonal to the images generated by 

the line vectors of \if{%) (30): 

n x (r) = i _ * (r ) (®t(r)#(r))^ *»(r) 



The estimate of T is given by (3 1 ) 

r = axgnunjjn^(r)Y|J 2 

The two alternative embodiments explained here above deal with signals 
sampled at output of the sensor or again for a signal correlated with a known signal 
having a value 1 . The matrices to be processed have a size NeNxP. 

According to another alternative implementation of the method which can be 
applied in the case of a calibrated antenna as well as that of a non-calibrated antenna, 
the method according to the invention uses correlated signals. 

Figure 4 shows a block diagram in which a correlator 10 receives, firstly, the 
signals received by the antennas and, secondly, the known signal. 

The correlator 10 has a number of outputs identical to the number of inputs, 
each signal si(t) that reaches an antenna being correlated by a known signal c(t). The 
different correlated signals si(t)*ci(t) where i is the index of the centre arc then 
transmitted to a device 1 1 adapted to the selection of a number of samples, especially 
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in order to work on a time interval defined beforehand. The signals that have been 
correlated and belong to the chosen time interval are then transmitted to a processing 
unit 12 comprising, for example, a microprocessor suited to estimating the values 9 
and % when the characteristics of the antenna are known or the value of % for an 
5 antenna having unknown characteristics. From these values, it is possible to 
determine the estimates of the impulse response of the propagation channel. 

The correlator 10 may be of the analog or digital type, the sampling steps 
being performed before or after the correlation step. 

The processing block also receives the sent signal, itself correlated for 
1G example by means of a correlator 1 3. 

When a known time interval is used, it is useful to know the time of reception 
at the block as well as the length of the channel. 

The signals considered in the method correspond to the signals output after 
correlation by the filter matched with the transmitted signals. 
1 5 Method for an antenna whose characteristics are known. 

In this part, the same assumptions are used as in the previous part explained 
in the case of an antenna for which the parameters are known. We shall again 
consider the signals to be known and shall examine only the case of the user o> in 
order to simplify the presentation* This model can immediately be extended to the 
20 MIMO (Multiple input Multiple Output) context. 

Modelling 

By correlating the signals received at each sensor with the signal of the user 
CO, we get (32) : 
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- Z *wk ft 1 ) - r v ) * ^H)) + w)) * n-t) 

By limiting the data processed to the part "containing" the propagation 
30 channel at output of the correlation, the complexity is greatly reduced because the 
size of the matrices used is diminished. 

By concatenating in the filtered samples at output of the sensors on a time 
interval Nl, the signal X can be modelled as follows (33) : 



35 
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where 



and 



x=.[y l5P (r e ),.-,y^(Mr e )J 
{(b(0 + b'(0) r *^'H)} jvr 



with *(tf,T)«[*i.--,*H 

and *,(«, *) - " *?> * ® a ^ 

The noise b to follows a normal law of covariance otRb where the matrix 
Rfo^R/'Ijy . R/* is the covariance matrix of the vector (34): 

The likelihood of the signal X 0 * is given by (35): 



P(X"mO, r), h«, 4, RJ = gCT 53 «* 



pvnv *r 



whence the log-likelihood (36): 



1 i 



L(X-|*(0,r),h- 4,£U) = -M^log(ir<4)-IoB(|Rb!)-^ |j R b s fX w - *(0,r)h«) 

By cancelling the derivatives of the log-likelihood with respect to a 2 and with 
respect to h w , we obtain the estimators of a* and h w (37). 

R ^(X"-$( g,r)hl" a 




and 



By replacing a 2 and h 10 by their expression in (36) we obtain the estimator sought tor 
the space-time parameters 9 and X (38); (39) 



0, t = arg min X w R£ l H£(ff, r)X 

0tT 



til 



where 



The estimated space-time parameters are used for example to determine the 
estimate of the impulse response of the propagation channel. 



Non-calibrated antenna 

In this part, we shall look at a reduced maximum likelihood method with an 
unknown antenna. The approach used will be identical to the preceding one. We 
shall give only the model of the signal, the log-likelihood and the estimators. 
5 Modelling 

The received signal is filtered by s m (-f) (40) : 



/UP* 

o n 



By concatenating the samples at output of the sensors and then the vectors 
obtained, we get (41); 

3 

Y<*=Scc* i +n to 



where 



y- = [if (t.), • • • ,< r (Mr e ),-.. ,xf(N,r.)] T 



with 
and 



The noise n w follows a standard law of co variance a^Fn where the matrix 
R„ — If, ® /? ?w . 

The likelihood of the signal Y* is given by (42): 
whence the log-likelihood (43) : 

L(Y w IS(r), a", 4, R n ) = -/V.iVlog^) - logflR,!) - 4" p£* (Y<* - S(r)cO 
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By cancelling the derivative of the log-likelihood with respect to <o\ and with 
respect to cc w , we obtain the estimators of a\ (43) and of a m (44): 



0 



md 5" - (S^R^))" 1 S(r)R^X" 

By replacing in (42) and by their expression, we find an 
estimator for the parameter T (45): 

r = axgmin Y^^H^rJY 1 - 

r 

where n s ^=I- S(r) (S(r)R^ »S(r)) ~ l 

The estimated value of the delay t is then used for example to estimate 
5 the value of the impulse response of the propagation channel. 

The different alternative modes of implementation of the method 
according to the invention, can be applied in many fields of signal processing. In 
communications systems, it is often necessary to know the propagation channel in 
order to demodulate the signals transmitted. 

D In view of spectral congestion and the growing demand from users, 

future systems will have several users communicating simultaneously in the same 
frequency band. Consequently, the techniques proposed here in the MIMO context 
have a direct application. 

Similarly, the technique of filter estimation with reduction of complexity 
can be applied in the MIMO and SIMO (Single input Multiple Output) context. 

These techniques are also applicable to radar processing. They can be 
used in a multiple-target or single-target context with echos 

Similarly, in radio altimetry, several echoes of the sent signal may be 
received and all of them have to be estimated. 
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Without departing from the framework of the invention, the invention 
may take account of the Doppler effect* For this purpose, one method consists for 
example in estimating the speeds by minimizing a multidimensional function that 
integrates all the parameters of the model, other than the complex amplitudes and the 
noise power, 

The method according to the invention can be applied especially in a 
UMTS telecommunications system comprising an array of antennas. 



